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Abstract—In this work, a robotic assistance strategy is 
developed to improve the safety in an artisanal task that involves 
a strong interaction between a machine-tool and an operator. 
Wood milling is chosen as a pilot task due to its importance in 
carpentry and its accidentogenic aspect. A physical model of the 
tooling process including a human is proposed and a simulator 
is thereafter developed to better understand situations that 
are dangerous for the craftsman. This simulator is validated 
with experiments on three subjects using an harmless mock- 
up. This validation shows the pertinence of the proposed control 
approach for the collaborative robot used to increase the safety 
of the task. 


I. INTRODUCTION 


A wood shaper is a machining tool used to profile wooden parts 
(see fig. 1). During the milling process, a craftsman is handling a 
non homogeneous anisotropic material, in strong interaction with 
the cutting tool. It can be difficult to guarantee the stability of the 
part throughout the task. This can lead to tragic accidents ranging 
from bruises to severe cuts of the carpenter’s limbs. According to 
a study conducted by J. Hester et al. [1], among 1200 apprentice 
carpenters interviewed, 80% of them have already observed a work 
accident. Similarly, in 2016, the French Institute for Craftsmen’s 
Health and Safety (IRIS-ST) counted 5118 accidents for 61764 
carpenters [2]. 27% of them were related to tooling tasks. 





Fig. 1: Illustration of a wood milling task by a trained expert 
on a wood shaper. The left arm maintains lateral stability 
while the right arm (hidden on this picture) pushes the part 
in direction of the cutting tool. 


Automated solutions for medium and large series production in 
the woodworking industry are not adapted to craft trades that re- 
quire flexibility and adaptability. A conventional robotic approach 
is too expensive for many small businesses and removes the human 
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from the loop. However, collaborative robotics potentially offers 
new possibilities for assisting carpenters, preserving their know- 
how and sensations while making the task safer. 

Several types of robotic assistance are generally envisioned [3]. 
One of the most commonly considered ones is comanipulation 
by gravity compensation and/or force amplification. For example, 
the authors of [4] propose that an upper limbs exoskeleton par- 
ticipates in reducing the strenuousness of some industrial tasks, 
the advantage of an exoskeleton being to act both on the posture 
while amplifying the user strength. This type of assistance is also 
considered in the construction industry for handling heavy loads 
[5]. In industrial domains where the workers are technical experts 
of complex manual tasks (i.e craftsmen), force amplification is 
also considered as a approach of interest. For example, Lamy et 
al. [6] have developed a robotic assistance solution for butcher 
cutting large pieces of meat in slaughter houses.While the impact 
on the worker is positive in terms of reduction of musculoskeletal 
disorders, the robot is not directly envisioned as way to protect the 
worker from potential cuts with the knife. 

While these research works provide examples of how collaborative 
robots can provide physical assistance to workers, this assistance 
is in general not designed with immediate safety in mind. In fact, 
prototyping robotic solutions dedicated to immediate safety in 
highly dangerous tasks is ethically complex since it can endanger 
subjects during tests. It then becomes difficult to analyze the tasks 
and design solutions without considering the use of simulation 
tools and/or meaningful but safe mock-ups. Simulation tools are 
well suited to systematically evaluate and optimize design choices 
related to the physical and control architecture of potential robotic 
assistant solutions[7]. 

In this work, a wood shaper simulation tool is proposed (sec- 
tion III). This simulator relies on a physical model of the system 
(section II) which accounts for the wrenches applied by the human 
operator, the cutting tool as well as a candidate robotic solution. 
This simulator is used especially in accidental situations. It is also 
used to evaluate the pertinence of a the proposed control solutions. 
The proposed approach to safety is to control the robot with the 
reduction of the amplitude of disturbances induced by a local 
modification of the nature of the wood (typically a node) in mind 
(section IV). A physical mock-up of the system is also proposed 
and the performances of the robotic solution obtained for different 
tuning of the controller are presented and analysed (section V). 


IT. TASK MODELING 


A. Equilibrium of the wooden part 


The work environment during a milling task consists of a cutting 
tool rotating at high speed around an axis, zpr, normal to an 
horizontal work table that serves as a support for the wood part 
and is defined by a plane P. A rail-shaped guide associated to an 
axis Tr € P is attached to the work table. 





Fig. 2: Schematic representation of a wood shaper. 


Let {R} = {O, TR, YR, ZR} represent a chosen orthonormal 
frame fixed in the workspace which origin is at the intersection of 
zr and plane P. 

The part is subject to a set of forces and moments applied 
by the operator and the cutting tool during the milling process. 
This corresponds to the more general situation of an object being 
subject to different external wrenches. The contact point associated 
(FT ml) (with f; € R3 and m; € 
R? the force and moment composing wrench w;) is denoted C; 
and is defined by the vector c; € R. A frame {C}; , with 
axes {£i, Yi, Zi} is also defined at this point. In the case of the 
wrenches applied by the operator, æ; is colinear to the human 
forearm axes and directed toward the object and y; is parallel to 
the support table plane. The nature of the wrench applied by the 
tool is introduced in section II-B. 

In order to derive the equilibrium conditions for the part, let 
express all wrenches at point O and in frame {R} 


to wrench w; = 


wiom=( Le, S ju (1) 
R;iS(Ci) Ri 
T 
Gi 


with S(C;) is the skew-symmetric matrix corresponding to the 
cross-product with vector OC; expressed in frame {C}, and R; 
the rotation matrix from frame {C}, to frame {R} . 

G; is the partial grasp matrix [8] which maps the contact wrench 
applied to the object at {C}, to the equivalent wrench expressed 
at O in frame {R}. The complete grasp matrix G € R°*°*"< 
of the system (where ne is the number of contacts) relate all 
external wrenches expressed at their respective contact points to 
the equivalent wrench expressed at O in frame {R} 


Wo/R = ( Gi Gn. )( wl = Wwy. Q) 
N 


G w 
The machining process results in a relative motion of the part with 
respect to frame {7}. In order to compute the motion resulting 
from the application of all external wrenches, one can apply the 


fundamental principle of dynamics to the system, using the grasp 
matrix defined in eq. (2) 


Gw = My + Bv. (3) 


where ~y and v are respectively the acceleration and twist of the 
part expressed at O in {R}. M € R®°*® and B € R°*® are 


respectively the symmetric and positive mass matrix and diagonal 
viscous friction matrix. 

Fore the sake of clarity and simplicity, and without loss of 
generality, it is assumed that the mass of the part is constant over 
the duration of the simulation, i.e the loss of mass due to the tooling 
process is neglected. It is also assumed that the friction between 
the table and the part is purely viscous. In practice dry friction is 
present but adding it to the model would not bring much added 
value to the overall proposed contribution. Finally, the weight of 
the part is fully compensated for by reaction forces of the table 
and any force applied to the part in the normal direction of the 
guiding rail is compensated for by reaction forces of the rail. 

As a consequence, in the unassisted case, the only external 
wrenches considered in the model are the ones of the operator’s 
left and right hand w; and wr and of the cutting tool we. 


B. Wood cutting model 


Several wood cutting models have been developed since the 
1950s. One of the first widely used models is proposed by E. Kivi- 
maa works [9]. However, this simplified model can lead to large 
differences (up to 40%) between the actual cutting forces and 
the modeled [10]. F Eyma et al. [11] have improved this model 
for 13 tropical wood species, but do not take into account the 
wood grain orientation. This characteristics largely impacts the 
magnitude of the cutting forces and, more importantly, a sudden 
modification of this orientation (typically with a node) represents a 
factor of potential accident. L. Cristovao et al. [12] account for this 
important effect in the wood cutting model. However their model 
only considers unidirectional cutting forces. 

The model used in this paper takes its origin in an experi- 
mentally validated model developed for medium and low density 
woods [10]. This model expresses the cutting force at the edge of 
the cutting tool and decomposes it in two components ftan and 
fnor which are respectively the cutting force tangential to the part 
feed velocity/trajectory (along the rail and aR) and normal to it 
(along yr). This yields the relation between the cutting wrench 
expressed at the edge of the cutting tool, the partial grasping matrix 
associated to the cutting tool and the cutting wrench expressed at 
O in frame {7} (index Le stands for cutting): 


We OJR = Gel Jtan Înor 01x4 ae (4) 


This model depends on several parameters including the cutting 
velocity and the density and orientation of the wood grain. These 
two last parameters are essential to simulate the appearance of a 
node. The expression of ftan and fnor taken from [10] are recalled 
here 


Jtan (a1+a2|cos(h;+a3)|"4)A1+(as+açs|sin(hs+a7)|"8)A2+a27 
fnor = (b1+b2|cos(dv+az)|°4) Bi +(bs+b6|sin(dv+b7)|°8)B2+ba7 
A: — a29 Ap? p21 y212 D713 Creer 1) 
Ag — 318 A p ° p720 ve21 D (asser 24) Hagge T 
Bı = ataa AO pb11yb12 p13 Cae e a 

Bo — a?18 719 9620 yb21 D2? (bag eo" 24) (byg e726) +, 


where ¢,: angle between the cutting velocity vector; ve and 
wood grains (rad) ; a, : thickness of the cutting layer also known 
as uncut chip thickness (mm) ; Ar : rake angle (rad) ; p : cutting 
edge dullness (um) ; ve : cutting velocity (m/s) ; D : wood density 
(kg/ m?) ; mc: moisture content (%); T : temperature of the wood 
(K) ; aj, by : coefficients determined by experimental identification 
in [10]. 


C. Modeling of the human operator 


Modeling the behaviour of human operator realizing an artisanal 
task such as wood milling is complex. In this work, several 
assumptions are made in order to obtain a model which level 
of complexity is appropriate given the exploratory nature of our 
study. While it does not seem desirable to include a complex 
musculoskeletal and motor control model of the human in this 
first study, the considered model should still account for the 
specificities of the task to be achieved as well as for some of the 
limits of the human in terms of applicable wrenches. 

Regarding the task specificities, expert teachers from a training 
centre for apprentice carpenters in Bordeaux, France (namely the 
BTP-CFA Gironde) have provided some general indications re- 
garding the use a wood shaper. First, the motion of the part during 
the cutting phase must be continuous, with constant velocity, in 
order to ensure a good cutting quality. Second the right hand 
(assumed to be the hand close to the guiding rail) is placed at 
the end of the part and provides wrenches mostly dedicated to 
its motion (along xR) while the left hand is placed on the side 
of the part and provides wrenches mostly dedicated to the lateral 
stabilization of the part (along yr) 

Regarding human capabilities in terms of producible wrench 
at the hands, the work in [13] provides information regarding 
the maximum force a male adult can produce with his hand in a 
pushing motion both in the direction of his forearm as well as in 
the direction tangential direction. The work in [14] also provides 
maximum values for the rate of variation of these forces induced 
by human motor reflexes and the dynamic of muscle contraction. 

Given these elements, the wrenches produced by the human 
operator during a wood milling task are computed in order to 
servo the velocity of the part in the feed direction aR given the 
other external wrenches applied to the part. Assuming that an 
experienced human operator trains a mental model of the part 
being milled and can estimate, through proprioception and contact 
with the part, all external wrenches at work, the ideal wrenches to 
produce using both hands can be computed using eq. (3) as: 


Gw, Wr + Gu, w] = M" + Bv — Gu.we (5) 


where: 

e Gw, and w; are the the partial grasping matrix for the wood 
shaper, the right an left hands (i = {c, r, l}) respectively; 

e v and wc are the twist of the part and the wrench applied by 
the cutting tool on the part respectively; 

e y is the control acceleration which is computed using a PI 
controller between the twist of the part v and the desired 
twist v*. The components of v* are 0 except in the Œ p linear 
direction which is the feed direction for the part. 

The bounds on the hand wrenches are written 


Wi min < wi < Wi mar», (6) 
Wi min < Wi < Wi max (7) 
where Wi min, Wi, maz, Wi,min and Wi max are taken from [13] 
and [14]. It is worth noticing that the hand pushing forces are uni- 


lateral (the operator can hardly pull the part) and the corresponding 
limit value f {7 r} min are set to 0. 


IHI. SIMULATION 
A. Simulation principles 


The aim of the simulator is to compute the evolution of the 
wrenches applied by the human operator on the part given a desired 


part velocity v*, the length of wood to be milled d for the part and 
a wood type. The initial state of the system is characterized by the 
initial positions of the hands in the reference frame {R} and by 
the initial twist of the part v(0) = 0. 

At each time step, the following computations are performed by 
the simulator: 


1) the partial grasp matrices (see eq. (1)) are updated given the 
current location of the hand contact points c; (i = {r,l}) 
with respect to frame {R}. Indeed these contact points 
remain constant in {C }; but are moving with respect to {R }; 

2) the tool cutting force is updated given the current twist of 
the part v; 

3) the PI controller is updated given the current velocity of the 
part v; 

4) the optimal hand wrenches are computed according to 
eqs. (5) to (7) (see hereafter); 

5) the part acceleration y is computed through eq. (3)) given 
the computed hand wrenches w?” (i = {r,1}), and tool 
cutting wrench we; 

6) the twist v of the part and location of the hands! c; (i = 
{r,l}) with respect to frame {R} are updated through 
numerical integration using an Euler integration scheme for 
the computed acceleration of the part +. 


B. Human operator wrenches optimization 


Computing the optimal hand wrenches (step 4) boils down to 
finding the optimal hand wrenches w?”* and w??" minimizing 
a cost function associated to eq. (5) under constraints defined by 
eqs. (6) and (7). 

Such an optimization problem can be written as a quadratic 
program under linear constraints (LQP) 


Wr,W] 
s.t. eq. (6) and eq. (7) 


where fw = ||Gw, wr+Gu wi -(M* + Bv — Gwu. we) || 
is the main cost function 

and freg = al ( w} wp E NW rece is the regularization 
cost. 

This part of the cost function is required for the optimization 
problem to have a unique solution. œ is chosen < 1. In practice, 
this leads to a solution both minimizing f as well as the weighted 
norm of the solution. The weighting matrix Wim is chosen to 
account for the bounds on the achievable human wrenches eq. (6) 
and favour directions where the force capabilities are the largest. 
Such a weighting matrix can be defined as a 12 x 12 diagonal 
matrix which elements are in R*: 











À 1 1 1 1 
diag = RE Gas . ©) 
tremax Mr, max Í nan MI max 


This simulator is implemented in Matlab and the LQP optimiza- 
tion problem is solved using the quadprog solver. The simulator 
handle 3D models but in the specific case of this planar application 
a 2D projection is performed. Fig. 3 provides a view of the visu- 
alization associated to the simulator. All parameters introduced in 
the paper can be tuned. However typical simulations are performed 
with the values provided in table I. 


lwhich is equivalent to updating the pose of the part. 





Fig. 3: Visualization of the force equilibrium between the 
hands and cutting forces for a desired velocity along £R of 
0.05m/s. 


TABLE I: Typical values of the simulation parameters (after 
2D projection) 


Simulator parameters 

Sampling time for integration 0.010s 

Maximum simulation duration 6s 

Maximum part travel distance 0.5m 

Desired feed rate/part velocity 0.05m/s 

Part characteristics 

Size (length, width, height) 

Mass matrix 

Viscous friction with the table 

Wood type 

Human operator characteristics 

Initial right hand contact point ( 

Initial left hand contact point ( 

Tool cutting edge coordinates ( 
( 
( 


(1.0,0.5,0.1) m 

diag(1.0kg, 1.0kg, 2.6042kg.m7) 
diag(30.0, 30.0,0.0) N.s.m7? 
Pinus Sylvestris [9] 


—0.95,0.20,0.0)2 m 
—0.80,0.45,0.0) m 
—0.05,0.05,0.0)2 m 

10.0, 10.0, 10.0) s—? 
10.0,10.0,10.0) s7? 

1076 

(160.0N, 15.0N, 13.0N.m)c,. 
(0.0N, —15.0N, —13.0N.m)c,. 
(13.0N, 133.0N, 13.0N.m)c, 
(—13.0N, 0.0N, —13.0N.m)c, 
120ms from Wi min t0 Wimax 


Proportional gains 

Integral gains 

Weight of the regularization cost 
Maximum right hand wrench 
Minimum right hand wrench 
Maximum left hand wrench 
Minimum right hand wrench 
Maximum wrench rate of change 


C. Simulating instabilities 


Two major cases of instabilities leading to accidents are gener- 
ally observed by experts. 

The first case of instability occurs when the cutting tool hits a 
node. Cutting forces suddenly increase and the carpenter’s limited 
reflexes and strength (eqs. (6) and (7)) may not be sufficient to 
instantaneously compensate for the unexpected increase in force. 
This leads to a transient inversion of the direction of the velocity 
of the part and an associated backward motion. The part hits the 
hands of the carpenter, potentially his arms and body. 

The second case of instability occurs when the cutting tool 
is in contact with a wood node and suddenly comes out of it. 
Cutting forces suddenly decrease and, due to similar limitations, 
the carpenter may not be able to instantaneously compensate for 
the unexpected increase in force. This leads to a transient increase 
of the velocity of the part and a faster motion of the part, and as a 
consequence of the hands of the carpenter, towards the cutting tool. 
The consequences can be dramatic, including complete severing of 
several fingers and/or deep cuts of the hands. Besides the more 
severe consequences on the human operator, this instability is 
harder to compensate for because of the unilateral nature of the 
constraint on the force for the right hand in the wz direction and 
of the very low force exertion capability of the left hand in the x; 


direction. 

Using this simulator, these two typical instabilities can be 
simulated by introducing nodes in the wood, i.e by varying the 
angle of the grain of the wood with respect to the feeding direction. 
Simulation results corresponding to such instabilities are presented 
in fig. 4 (a) and (b) for the first case of instability (cutting tool 
entering a node) and fig. 4 (c) and (d) for the second one (cutting 
tool leaving a node). 

On fig. 4 (a) and (b), slightly before t = 3.1s, the simulated 
apparition of a node leads, as expected, to an instantaneous” 
increase of the cutting force from ftan = —39.8N to ftan = 
—132.2N. This induces an increase of the required force f; for the 
right hand in the direction of motion of the part (along ær). The 
limit on the rate of evolution of the wrench applicable by the hands 
is such that the part experiences a fast backward motion ( 25cm 
and peak at v” < —1.0m/s) before getting back on track once 
the hands have reached the required level of force. Similarly, on 
fig. 4 (c) and (d), the simulated disappearance of a node leads to an 
instantaneous decrease of the cutting force from fran = —112.0N 
to ftan = —20N. Analogously to the previous case, this generates 
a fast forward motion (same amplitude) of the part before the 
motion of the part is performed at constant velocity again. 

These two examples confirm the ability of the simulator to 
reproduce physically realist simulations of both nominal and ac- 
cidental situations. 


IV. ROBOTIC ASSISTANCE FOR SAFETY 


The robotic assistance provided to the user should aim at 
decreasing the risk of sudden modification of the velocity of the 
part while allowing to maintain direct physical interaction of the 
carpenter with the wood as well as his sensations. Candidate 
control modes should also imply the least possible need for exter- 
nal sensing devices which may not be robust enough in a dusty 
environment and largely increase the complexity of the overall 
solution. 

To meet this specifications, it is proposed to use the robotic 
system in order to damp the transient velocity change when 
an instability occurs. The robot end-effector is assumed to be 
maintained in contact with the wood (the study of the associated 
technological solutions is out of the scope of this work) at an 
arbitrarily chosen point C'a (index a la stands for assistance). 
The associated wrench wa and associated partial grasping matrix 
Ga can be included seamlessly in eqs. (3) and (5) and section III-B 
and the simulation directly includes the contribution of the robot 
to the total wrench applied to the part. 

Assuming the use of a torque control robot (e.g. 7DOF Panda 
from Franka Emika), the emulation of a damping behaviour can 
be obtained by choosing the control torques T° as 


T° = g(q) — J(q) Bat (10) 


where, for the chosen robot: t is the twist of the end-effector, 
directly related to v, the twist of the part ; q is the generalized 
coordinates vector ; g(q) is the gravity compensation torque ; J(q) 
is the end-effector Jacobian matrix and B, is a damping matrix. 
This matrix is typically diagonal with components defined in R*). 
The damping coefficients can for example be chosen to be higher 
in Cartesian directions requiring more damping. In this work they 
are chosen, without loss of generality, to be equal (ba). In the 


*in practice this increase is not instantaneous but can be considered as 
such given the reaction time of the human. 
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Fig. 4: Results of the simulation of an instability due to the 
cutting tool entering ((a) and (b)) and leaving ((c) and (d)) a 
node. Figures (a) and (c) plot the evolution with respect to 
time of the distance travelled by the part along a as well as 
the corresponding part velocity in the same direction. Figures 
(b) and (d) provide the evolution of the right and left hands 
forces as well as of the cutting forces in the XY plane of 
the reference frame. 


context of this work these directions are obviously xR and yr. 
With this controller, the robot wrench applied to the part, is given 
in the robot end-effector frame by 


V. EXPERIMENTATION AND RESULTS 


To confirm the observations made in simulation on the impact 
of a sudden variation of the cutting force on the part motion, an 
experiment has been designed to reproduce one of the two cases 
of instability described in section HI-C. It was considered more 
interesting to focus on the case where the cutting tool exits a 
node given the higher risk of induced dramatic accident (potential 
contact of the hands with the cutting tool). 

While conducting experiments using a real wood shaper would 
have been possible, it was judged safer to realize this first experi- 


ment using a simplified mock-up that cannot cause any harm to the 
user (see fig. 5). 


A. Mock up description 


Fig. 5: Experimental mock-up including a table, a guiding 
rail, a wooden board, the Panda robot and one of the three 
subjects. The pulley system and load used to produce the 
emulated cutting force are placed under the table. 


To remove any source of risk, the cutting force is emulated using 
a constant horizontal force of 75 N applied to a wooden board. This 
force is produced using a load attached via a pulley to the board. 
So, for the wooden board to move forward, the user must apply a 
pushing force of at least 75N. The motion of the board is guided 
by arail similar to the ones used on real wood shapers. The sudden 
force drop due to the cutting tool exiting the node is emulated by 
an instantaneous reduction of the dragging force to 25N. Note that 
to guarantee the safety of the test subjects, the force implied in the 
process are reduced compared to the more realistic cases simulated 
and presented in section II-C. 

To verify the hypothesised benefits of a robotic assistance, a 
7-DOF torque controlled robot (Panda from Franka Emika) is 
attached to the wooden part. The control law described by eq. (10) 
is used at all time with damping values in the XY plane of the table 
ranging from 0 to 600N s.m |. The robot being rigidly attached 
to the wooden board, it provides access to its position and velocity. 


B. Protocol 


Three healthy male subjects (none of them is co-author of this 
paper) are asked to move the wooden board part along the rail at 
a constant velocity. Indications on where to position their hands 
are given to them. The subjects first perform a few trials without 
accident emulation to get accustomed to the dynamics of the 
wooden board subjected to the emulated cutting force. After this 
short training, each subject realizes five recorded experiments with 
accident emulation, each time with a different value of the damp- 
ing coefficient. For subject 3, an additional downwards constant 
vertical force of 50N is applied by the robot to the wood when the 
rate of change of the planar velocity is greater than 20m/ s?. 


C. Results and discussions 


Figure 6 provides several information regarding the performed 
experiments. Figure 6 (a) (blue curve) compares the velocity 
profile of the part in the £r direction when the disturbance occurs 
for subject 1 for different damping values. The effect of damping 
on the reduction of the peak velocity is clearly positive in terms of 
reduction of the forward motion of the hands induced by the wood 
disturbance. Interestingly, considering a similar type of instability 
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Fig. 6: (a) shows the velocity evolution of the wooden part for different damping values provided by the robot (from the 
unassisted case ba = ON :s-m 1 to the maximum damping applied value ba = 600N -s-m +). The blue curve represents the 
data measured during the experiment of the first subject and the red curve the one generated by the simulator. (b) describes 
the evolution of the velocity of the wooden part for subject 3, i.e. in the case where a 50N vertical force is also triggered 


in case of disturbance (red and green curves). The blue curve corresponds to the fully unassisted case. 


in simulation leads to a similar effect (red curve), both qualitatively 
and quantitatively. As expected, the larger the damping value, the 
more efficient the reduction of the peak velocity is. Nevertheless, 
one should note that, quite logically, increasing the damping has 
a slightly noticeable effect on the effort required by the simu- 
lated human to maintain the part velocity at the desired value: 
the steady-state force to be applied by the right hand along ap 
increases from 72 to 86N when the damping increases from 0 to 
600N.s.m—*. 

Finally, fig. 6 (b) presents some of the results for subject 3 
for which the control action of the robot is both a damper in the 
planar direction and a brake in the vertical direction. This constant 
braking action which is tuned to be triggered in case of abrupt 
acceleration induce in practice an increase of the friction between 
the wooden board and the table and thus reinforces the primary 
damping effect. Three cases are exposed: the completely unas- 
sisted one (blue curve), the one where only the vertical breaking 
force is applied (red curve) and the one where maximum damping 
is applied in the planar direction additionally to the vertical braking 
action. While the case of maximum planar damping is rather 
comparable to the results obtained with subject 1, the pure vertical 
braking action (red curve) produces, quite interestingly, a very 
efficient reduction of the effect of the wood disturbance. 

The overall obtained results are summarized in table II. This 
results are in-accordance with the conclusions drawn for fig. 6 and 
also confirm that the addition of a vertical force significantly re- 
duces disturbance peaks. This is of interest as this vertical braking 
action 1) may not require the robot to be attached to the part (this 
is clearly an advantage in terms of applicability) and 2) is only 
triggered in case of disturbance and does not induce a permanent 
increase of the required pushing force of the operator. It may 
however be more complex to tune than the damping action. Further 
comparisons of these two strategies require a more longitudinal 
study which could be the topic of some future work. 


VI. CONCLUSION 


This article presents a robotic assistance solution that secures a 
wood milling task, well known as dangerous. A simulated physical 
model shows the interaction wrenches between the human, the tool 
and the wooden part as well as their dynamics during the task. 
This model allows to analyse accidentogenic cases and is also 
used to design a proper control strategy for the assistive robot. 
The proposed approach is validated experimentally by testing, on 


TABLE II: Effect of the damping coefficient on the velocity 
peak for a desired velocity of the part v* -£r = 0.025m/s. 
The velocity peak is expressed in m/s. The relative reduction 
(%) of the peak height with respect to the the unassisted case 
(ba = ON.s.m t) is also provided. The results are presented 
both for the simulation and for the three subjects who per- 
formed the experiment. An additional 50N vertical braking 
action is triggered for Subject 3 when the acceleration of the 
part is greater than 20m/s°. For this subject, the value of 
the velocity peak in the case with no damping is provided 
for the fully unassisted case, 1.e. the vertical force is not 
triggered. The value of the velocity peak corresponding to 
zero damping and the triggered vertical force is 0.15m/s 
which yields a reduction of 42.3% of the peak. 


ba | 0 [100 % [200 % |400 % [600 % 


Sim. 0.10 67.71% 
SI 0.10 65.5% 
52 0.12 53.8% 
53 0.07 73.1% 





a 7-DOF robot, a Cartesian space damping controller to secure 
an emulated accident on an harmless mock-up for three healthy 
subjects. The experimentation confirms the results obtained in 
simulation and the pertinence of both the simulator and the control 
approach. 

These promising results suggest that simulation can be used 
to design advanced control architectures including virtual walls 
around the cutting-tool as well as adaptive damping that would 
account for the operator know-how and level of expertise. Im- 
provement of the simulated human model could also permit the 
design of the assistance solution with ergonomics in mind. 

Overall, this work brings to light the artisanal field that can 
clearly gain a lot in terms of safety and worker’s health preser- 
vation through the use collaborative robots. 
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